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The electronic spectrum on the spherical surface of a topological insulator reflects an active 
property of the helical surface state that stems from a constraint on its spin on a curved surface. 
The induced spin connection can be interpreted as an effective vector potential associated with 
a fictitious magnetic monopole induced at the center of the sphere. The strength of the induced 
magnetic monopole is found to be gr = ±27r, being the smallest finite (absolute) value compatible 
with the Dirac quantization condition. We have established an explicit correspondence between 
the bulk Hamiltonian and the effective Dirac operator on the curved spherical surface. An explicit 
construction of the surface spinor wave functions implies a rich spin texture possibly realized on 
the surface of topological insulator nanoparticles. The electronic spectrum inferred by the obtained 
effective surface Dirac theory, confirmed also by the bulk tight-binding calculation, suggests a specific 
photo absorption/emission spectrum of such nanoparticles. 



I. INTRODUCTION 

It was only several years ago that the idea of topolog- 
ical insulator has been proposed as a possible candidate 
for the new state of matter in the field of condensed- 
matter ii The original theoretical idea has already been 
extended in various aspects, made applicable to a broader 
range of phenomena, including superconductivity and 
superfluidity.— The related research areas arc now re- 
classified and recognized as that of the topological quan- 
tum phenomena. Naturally, the outbreak of this new re- 
search field owes much to a rapid success of experimental 
studies that have demonstrated that the new theoretical 
idea has much reality.— 

The existence of a single gapless Dirac cone in its sur- 
face spectrum is a hallmark of strong topological insu- 
lators. Here, we focus on a specific property of this ro- 
bust and protected surface state on a curved surface,— 
the " spin-to-surface locking" . It is indeed specific to the 
topological insulator surface state and distinguish it from 
other realizations of gapless Dirac cones in condensed 
matter such as in graphene^i^ and related carbon mate- 
rials. The role of spin-to-surface locking may be most 
accentuated in the (pseudo-cylindrical) wire-shaped ge- 
ometry in which an anomalous Aharonov-Bohm type of 
oscillation has been reported^ Motivated by the reality 
of such transport measurements which may allow for a 
direct observation of the spin Berry phase, theorists have 
extensively studied the role of this phenomenon in the 
transport characteristics of the surface state.—"— 

A remarkable consequence of the spin-to-surface lock- 
ing in the cylindrical geometry is the half-integral quan- 
tization of the orbital angular momentum. Clearly, such 
half-integral quantization leads to appearance of a finite- 
size energy gap in the hitherto gapless surface electronic 
spectrum. Interestingly, introduction of a physical mag- 
netic fiux of half of a unit fiux quantum through (pierc- 
ing) the cylinder compensates the Berry phase associ- 
ated with the spin-to-surface locking, and closes the gap. 



The same mechanism applies to the classification of gap- 
less electronic states bound to a crystal dislocation line 
penetrating an otherwise surfaceless sample of a three- 
dimensional topological insulator.— A more systematic 
consideratioE^i^ on such gapless electronic states associ- 
ated with a topological defect in a topological mother 
system has been developed from the viewpoint of clas- 
sifying topological insulators and superconductors in a 
unified way solely from their symmetry classJ^"— 

The specificity of the cylindrical surface is that it is flat 
in the sense that it has everywhere a vanishing Gaussian 
curvature. On the surface of a topological insulator of 
more generic shape or geometry yielding a finite curva- 
ture, the effect of spin-to-surface locking mentioned ear- 
lier will be modified by that of a finite curvature. A 
spherical surface of topological insulator— is a prototyp- 
ical example in which such an interplay is expected. We 
show in this paper that the two effects are both expressed 
in terms of a Berry phase, but of contrasting nature (see 
Table |l| . The two types of Berry phase both contribute 
to the formation of a finite-size energy gap. The result- 
ing surface electronic spectrum on the sphere is shown to 
have a substantial compatibility with the result of tight- 
binding calculations performed for a cubic system (for 
tight-binding calculation involving the bulk, cubic im- 
plementation is much straightforward). A related but 
different scenario on the fate of such a (planar) gapless 
Dirac cone embedded on the curved spherical surface has 
been proposed in the study of the electronic states in 
fullerenei^"— 

In addition to the spectrum, the structure of the sur- 
face spinor wave function is another highlight of the pa- 
per. On the curved spherical surface of a topological insu- 
lator the strong spin-orbit coupling in the bulk, combined 
with the twisting of the phase shift due to the two types 
of Berry curvature, leads to a non-trivial spin texture. By 
explicitly constructing the surface spinor wave function 
we reveal such a rich spin texture possibly realized on the 
surface of topological insulator nano-particles. 

The paper is organized as follows. In Sec. II, the 
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effective surface Dirac tlieory is derived from tlie gapped 
bulk Hamiltonian, in which two types of Berry phase 
appear. The nature of these two types of Berry phase is 
discussed and contrasted in Sec. III. The solution of the 
effective surface Dirac equation is given explicitly in Sec. 
IV. The surface wave function is shown to be expressed in 
terms of the Jacobi's polynomials. The obtained discrete 
energy spectrum is compared with the result of (bulk) 
tight-binding calculation in Sec. V. This leads us to our 
conclusions. Some details of the formulation are left to 
the appendices. 



II. DERIVATION OF THE SURFACE 
EFFECTIVE HAMILTONIAN 

Let us first derive an effective "Dirac operator" on 
the spherical surface, starting with a bulk Hamiltonian. 
Our starting point is the following gapped bulk effective 
Hamiltoniar^i^l in the continuum limit: 



bulk 



e{p)l + ni{p)T^+ AT^{p^(Jx+Py(7y+p^(j^), ( 1 ) 



describing a three-dimensional (3D) Z2 topological insu- 
lator, where 



m{p) = mo + m2{pl +pI+pI), 



(2) 



is a (generalized) mass term containing both the constant 
and quadratic (Wilson) terms. For simplicity, we have 
chosen the Wilson term to be isotropic. The two types of 
Pauli matrices cr ~ (o"a;, Cy, CTz) and r = {Tx,Ty,Tz) rep- 
resent, respectively, the real and orbital spin degrees of 
freedom, and 1 is the 4x4 identity matrix. The Hamil- 
tonian ([T]) is time-reversal invariant, i.e., invariant under 
the operation of = iayK, where K represents complex 
conjugation. The two types of Pauli matrices represent 
independent degrees of freedom acting on spinors living 
in a different space. To make this point explicit one may 
express Eq. ([1]) in the following 4x4 matricial form. 



i^buik = e(p)l + 



m(p) Apz Ap- 

Apz —m{p) Ap- 

Ap+ ■m{p) -Ap:, 

Ap+ —Ap, —m{p) 



, (3) 



where p± = Px i ipy ■ 

In the following demonstration, we choose e{p) to be 
null so that the spectrum be symmetric with respect to 
E = (particle-hole symmteric). The vanishing of the 
e(p)l term upgrades the symmetry of the model from 
class AH to Dili (see Appendix A for details), but leaves 
unchanged the distinction between topologically trivial 
{mo/ 1712 > 0) and non-trivial (too/™2 > 0) phases. Note 
that the minimal model we consider contains only three 
control parameters, mg, fTi2 and A. We also consider the 
spherical geometry, assuming that the topological insula- 
tor described by Eqs. ([T]) and ([2]) occupies the interior of 
a sphere of radius R. We introduce standard 3D spherical 



coordinates: (r, 6, 
as 



related to the cartesian coordinates 



X = r sin cos < 
y = rsinO sin q 
z = r cos 0. 



(4) 



The momentum components, {px,Py,Pz) hi Eq. ([3]) ex- 
pressed in the cartesian coordinates are rewritten in 
terms of the derivatives with respect to the spherical co- 
ordinates, (r, 0, 0), by following the standard procedure. 
We also introduce unit vectors, r, 6 and <p, pointed, re- 
spectively, in the direction of the increase of (r, 9, 4>). The 
momentum operator p = p^x + pyij + p^z can be re- 
projected onto the directions of such unit vectors in the 
spherical coordinates as p = p^r + pg9 -\- Ptpfj), where 
Pr = -i{dr + 1/r), pe = -ide/r and p^ = ~id^/{r sinB). 

To derive the surface effective Hamiltonian in the spirit 
of k ■ p-approximation, we divide i/buik h^to two parts: 
^^bulk = H± + i^ii , _ 
solve the radial eigenvlaue problem 



where H± = iJ|pa=o.p^=Oi and first 



(5) 



instead of -ffbuikl^) = E\ip). Let us consider a possible 
form of the surface solutions of Eq. ([S])- IV') ™.ay take 
the following form: 



\^P) = \tP{r,e,cl,)) 



■R)\ 



t(0, </>)), 



(6) 



where k"^ measures the penetration of the surface wave 
function into the bulk. Here, taking a linear combination 
of the solutions of the form of Eq. ^ we construct a so- 
lution of Eq. ([5]) which is compatible with the boundary 
conditionJ^i2^"i2^ 



\'<pir = R)) = 0, 



(7) 



i.e., the condition that all the four components of the 
wave function -0 vanish on the surface of the sphere (at 
r = R). As shown in Appendix B, this can be matched 
by superposing two damped solutions of the form of Eq. 

Importantly, the solutions of such a boundary value 
problem must satisfy the zero-energy condition (jB14p . 
i.e., the (surface) Dirac point is at £' = in our model. 

The zero-energy condition (|B14I) helps simplify the 
solution of the radial eigenvalue problem ([5]). Since 
E±_ — 0, solving Eq. ([5]) is equivalent to finding that 
satisfiesr^i^ 

r.H±|7/0=O, (8) 

where for taking the form of Eq. ([6|), 

TzHi^\4,) = [m^ [k] + TyAnf • cr] | V) , (9) 

where ?7!,^[k] = mg — TO2K^. This implies that the or- 
bital part of the eigenspinor u(0, <f)) can be chosen as an 
eigenstate of - 



T-ykyi) = ±\Ty±). 



(10) 
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TABLE I. Characterization of the two types of Berry phase manifesting on the surface of a spherical topological insulator 



type 


(A) 


(B) 


(geometrical) origin 


curvature in the polar (9-) direction; effect of 


effect of rolling the surface in the azimuthal 




closing the surface at the north and south poles 


{(j)-) direction 


appearance 


in the covariant derivatives, or 


—id^ — ^ —id^ + 1/2 


(where, how) 


de^de + ^ cot (61/2) 




shifting the spectrum? 


yes 


yes 


relation to spin-to-surface breaks the locking 


expression of the (tendency to) 


locking 




spin-to-surface locking 


sensitivity to the choice 


no 


yes 


of basis 






other examples? 


fuUerene (buckyball) 


cylindrical TI 



To fix the notation, let us express the explicit vectorial 
representation of |tj,±) as 



1 

71 



1 

71 



1 



(11) 



The real spin part of Eq. (jH]) can be also diagonalized by 
pointing the eigenstates of cr^ in the direction of f, i.e., 
by 



r- = 



1 

71 
1 

71 



cos ■ 



i4>/2 



am ■ 
cos ■ 



(12) 



Combining these two types of spinors, one can compose 
the spinorial part of ip that can be matched with the 
condition ([5]), leading to 



[mj_[K] + KAa/3]\ra)i3 = 0, 



where 



Eq. ([T5)) implies 



-al3A ± + 4too"12 
2m2 



(13) 



(14) 



(15) 



The radial part of the wave function p{r) that is compati- 
ble with the boundary condition ([7]) takes the form given 
in Eq. ()B9|) . here, with k± being the two solutions of 
Eq. (jlSp . Clearly, both and k_ must be positive, for 
the wave function ip to describe a surface state localized 
in the vicinity of the spherical surface (r ~ R). Thus, 
in order to cope with the boundary condition, one must 
have both a/3 < and 11107712 < for the choice of model 
parameters such that A/m2 > 0. Notice that the second 



condition, mom2 < 0, which has appeared here automat- 
ically from the boundary condition, is a requirement for 
the system to be in the topologically non-trivial phase. 

We have thus successfully found the two basis eigen- 
states of H± for constructing the effective surface Hamil- 
tonian. For simplicity of the notation we denote them as 
-0 = |±)), where 



-)) =p(r)|f+)_, 
-))=p(r)|f-) + . 



(16) 



To avoid misunderstanding of the notations let us express 
explicitly the four-component vectorial form of |r±)zp 



1 

—i 
1 



e">'/'^ COS ^ 



(17) 



Here, the arrangement of the basis is made in accordance 
with that of Eq. Notice that the eigenvectors of Eq. 
(|17p are double-valued with respect to the azimuthal angle 
(j). This does not happen for the polar angle 6, since the 
domain of definition for 9 is restricted to a finite range, 
9 €z [0, tt], and not periodic, in contrast to (j). The double- 
valuedness stems from our choice of the (arbitrary) phase 
factor in front of Eq. (|T2)) . This is, on the other hand, 
merely a choice, and one can equally formulate the same 
problem consistently using a pair of single-valued eigen- 
vectors. We leave further arguments on this point to Sec. 
Ill and here take these double-valued eigenvectors as a 
basis for constructing the surface effective Hamiltonian. 
The effective surface " Hamiltonian" Hdv acts on a two- 
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component spinor, 



(18) 



Within the fc -p-approximation any surface state \a)) can 
be represented as a hnear combination of |+)) and |— )) 
with the amphtude specified, respectively, by a+ and q;_, 



Apart from an overah constant in front of the expression, 
this can be identified as the "Dirac operator" for a free 
massless fermion on the sphere.—"—'^"— 

The origin of the Berry phase term can be attributed to 
the covariance of the derivatives dg and 9^ on a curved 
spherical surfacei^"— In this regard, the Berry phase 
term appears as a spin connection in I?^ as 



|a))=a+|+))+a_|-)) 



(19) 



and 



= dij, + cos 9, 



replacing in Eq. (f24| as 



(25) 



— 1-LdvOL- 



(20) 



The explicit form of T-L^v can be determined by evalu- 
ating each of the matrix elements H\\ against the basis 
vectors |±)), i.e., ((±|_ff|| |=f)). The procedure we follow 
here is precisely in parallel with that of the standard 
degenerate perturbation theory. = H± is an unper- 
turbed Hamiltonian and |±)) arc its (two-fold) degener- 
ate eigenstates. To find the (degeneracy-lifted) spectrum 
of the perturbed Hamiltonian, iftot = -^^o + H' , where 
H' = H\\ , i/tot = ^^buik we first calculate the matrix 
elements: {{a\H'\fi)) = {Hdv)a[i («,/? = ±), then diago- 
nalize this 2x2 coefficient matrix. 
The explicit matrix form of H\\ is 



m|| —lADz 
'lADz ~™|| 
-iAD+ 



-iAD 



+ 







-iAD^ 
-iAD^ 
m|| iADz 

— m|| 



iAD,, 



, (21) 



where D± and are defined in Eqs. (|Bip and (|B2p . 
Performing the r-integral in ((a|7J|| |/3)) one can safely 
replace the r-dependence in these expressions with the 
radius R of the sphere, assuming that the surface wave 
function is well localized in the vicinity of the surface. 
Alternatively, one can equally regard D± and as 



R 



d 

cose-± 



d 



smf 



Dz 



d 



R do 



(22) 

At leading order in the expansion with respect to 1/i? 
the diagonal terms of Hn can be neglected, since 



mil 



m2- 



1 

i?2 



(23) 



[see also Eqs. (|B3P and (jB5[) ]. Within this accuracy the 
coefficient matrix "Hdv is found, after some algebra, to be 





sin 6 



-dg 



cot g 

2 



cot t 
2 



(24) 



-idg 
R 



RsinO 



(26) 



Alternatively, the Berry phase term can be absorbed in 
dg by introducing 



dg ^dg + - cot ( 



(27) 



In terms of dg , the Dirac operator ((24)) can be also rewrit- 
ten as (c/. Table I) 



Hdv--^ 

^A 



-dg 
sine 



sin I 



-idg 
R 



-id, 



'0 



i?sin( 



(28) 



III. 



NATURE OF THE TWO TYPES OF BERRY 
PHASE 



The advantage of considering the spherical geometry 
is that the existence of two different types of Berry phase 
becomes apparent; each associated, respectively, with an 
electonic motion in the polar (6-) [type (A)] and az- 
imuthal (0-) [type (B)] directions (see Table|T|. The type 
(A) Berry phase is intrinsic to the curvature of the spher- 
ical surface, while the type (B) is associated with the so- 
called spin-to-surface locking . On a cylindrical surface, 
on contrary, only the latter [type (B)] manifests, since 
the cylindrical surface has a vanishing Gaussian curva- 
ture. The contrasting behaviors of the two types Berry 
phase are summarized in Table H] 

To highlight the distinct behaviors of the two types of 
Berry phase, let us reconsider the Dirac operator (p4)) 
expressed against the double- valued basis vectors PT| . 
As mentioned earlier, this was not a unique choice of the 
basis. One can equally choose them to be single-valued 
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as 



vector potential, e.g., 



1 

—i 
1 

—i 

1 

i 



cos : 



e ^ sm ■ 



sm ■ 



e*"^ cos I 



(29) 



This type of a single- valued choice of the basis is often 
employed in the fc-p-description of the electronic states in 
graphencii Once this choice of basis vectors is adopted, 
one can repeat the same procedure as we described in 
the last section, to find the surface effective Hamiltonian 
"Hsv, or the Dirac operator in this basis, as 



A 
R 



2 2 



_ 1 cot ^ 
sine 2 2 





(30) 

In passing from Eq. ((24)) to ([30]) . the matrix elements are 
replaced as 



(31) 



Here, the additive factor 1/2 is nothing but the Berry 
phase TT associated with the spin-to-surface locking in 
cylindrical surfaces that has been extensively discussed 
in the literature<^"— 

The appearance of the type (A) Berry phase (see Ta- 
ble m is not restricted to the topological insulator sur- 
face state. It has already appeared in the study of 
the electronic spectrum of fuUerene, typically the one 
called, "buckyball" (or Buckminsterfullerene) <^"— The 
type (B) Berry phase, on contrary, is specific to the topo- 
logical insulator surface state. By its nature whether this 
type of Berry phase appears explicitly in the effective 
Dirac operator depends on the choice of the basis. In 
Hdv [Eq. the Berry phase is superficially hidden in 

the anti-periodicity of the basis spinor (|17|) with respect 
to the azimuthal angle (j). When one considers the or- 
bital part of the wave function or spinor (|18|) . this point 
must be carefully taken into account in its periodicity 
with respect to (j). This point will be clarified in the next 
section. 

The Berry phase term, or more precisely the spin con- 
nection of the form of ± cot 0/2 in Eq. ([24]) . or equiv- 
alently, either cot(e'/2)/2 or tan(6'/2)/2 in the two off- 
diagonals of Eq. (pO| . can be interpreted as a vector 
potential generated by an effective magnetic monopole. 
Indeed, the magnetic field associated with a magnetic 
monopole of strength g can be successfully encoded in a 



Ai = - — tan - cj), 

Anr 2 

-9 S 1 
Ajj = — — cot - d>, 



(32) 
(33) 



by introducing the concept of Dirac's string. Here, Aj 
and Ajj correspond to a choice of the gauge in which 
the Dirac's string runs, respectively, on the — z- {+z-) 
axis [direction of the south {9 = tt) vs. north {9 = 0) 
poles]. Eqs. dH]), imply that the strength of the 
"induced" monopole is, respectively, g = —2tt for 
and g = 2tt for a_. A fictitious magnetic monopole of an 
opposite charge is effectively induced at the center of the 
sphere for the two spin components of the surface spinor 
wave function [see also Eq. (01])]. 



IV. SURFACE EIGENSTATES ON THE SPHERE 

To find the electronic spectrum of the spherical topo- 
logical insulator surface state, one needs to solve the sur- 
face Dirac equation. 



(34) 



explicitly. In accordance with the decomposition of the 
Hamiltonian: i?buik = H± + , we have also separated 
the energy eigenvalue E of the original eigenvalue equa- 
tion. 



-f^buikla)) = E\oi)) 



(35) 



into two parts: E = E± + E\\. But, of course, since E± = 
[Eq. (|B14[) ]. E = E^^. As we have already emphasized, 
there exists some freedom in the choice of basis vectors, 
but here mainly to ease comparison with the literature 
we take them double- valued as in Eq. (fTT]) . To solve 
the eigenvalue equation (|34p explicitly we first separate 
the variables as a.{9, (j>) = e^""f'a.m{9), or in terms of the 
components, 



(36) 



Clearly, the quantum number m signifies physically the 
z-componcnt of the orbital angular momentum. Here, 
the important point is that this m takes half-odd integral 
values, 



a+(0,0) 






_a-(^»_ 




a„i-{9) 



2' 2 2' 



(37) 



in contrast to the case of fullerene i^^i^^'^^ Any solution of 
the original Scrodinger equation -ffbuikjV') = E!\ip), or its 
surface solution \a.)) [sec Eq. ([TO]) ] obeys to a periodic 
boundary condition with respect to (j). To be explicit, 
any bulk (or surface) solution \ip) = \ip{r,9,(j))) satisfies 
Itpir, 9,4> + 27r)) = \ip{r, 9, (f>)), and in the same sense. 



\ai9,4> + 27r))) ^\a{9,^))). 



(38) 
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On the other hand, the basis spinor (jl7|) spanning the 
space of the surface solutions | a (6*, </>))) [Eq. p^ ] is 
anti-periodic with respect to 27r-rotation of the azimuthal 
angle, i.e., |±)) — under (j) ^ (f> + 2tt. In or- 
der to insure the periodicity of the total wave func- 
tion p8|l this minus sign must be compensated, or ab- 
sorbed in the prefactor a± = a±{9,(t)) of Eq. ((T^ . i.e., 
a±(6',0-H27r) = -a±(6l,0) 

In terms of a„i+(6') and a„i_(0), the eigenvalue equa- 
tion (p4)) becomes a couple of first-order linear differential 
equation, 



A 


■ d 


m 


cote 


R 


d0 ^ 


smO 


+ 2 


A 


■ d 


m 


cote 


R 


d9 ^ 


sin 6 


+ 2 



These two equations combine to give 
Id d 1 



smede^^"^ de 



+ x 



sm 

2 



m cos ( 

2 



(40) 



in which two important parameters a and A have been in- 
troduced. CT = ± specifies the spin index in the subscript 
of amcr, whereas A parametrizes the energy as 



En 



A. 
R' 



(41) 



Let us remark here that Eq. (|40p is formally equiv- 
alent to a differential equation defining the so-called 
monopole harmonicsj^l the latter describing (the angular 
part of) the electronic motion in the presence of a mag- 
netic monopole. The role of the non-trivial Berry cur- 
vature imposed by the spherical geometry can be thus 
interpreted as an effective magnetic monopole induced at 
the center of the sphere. An electron in the surface state 
of a spherical topological insulator "sees" such an effec- 
tive magnetic monopole and behaves accordingly. In the 
notation of Ref?^, the quantum number q specifics phys- 
ically the strength g of the effective magnetic monopole 
placed at the origin as g = A'Kq. The Dirac's quantization 
condition restricts the "allowed" value of this quantum 
number q to be half integral: <? = 0, ±1/2, ±1, • • • . Here, 
in Eq. pO|) . q is identified as 



2 



±- 



(42) 



The cigcnfunction amaiO) of Eq. (|40| is indeed related 
to the monopole harmonics specified by this value of q. 
Thus at the center of a spherical topological insulator 
a magnetic monopole of strength ±27r, the smallest fi- 
nite value compatible with the Dirac's quantization con- 
dition is effectively induced. By its nature the induced 
monopole is automatically regularized by a Dirac's string. 



Introducing a new independent variable C = cos 0, one 
can rewrite Eq. as^^ 



dC dC 



- am( + 1/4 



1 



+ A^ - - ama = 0. 

This can be further rewritten in the form of Jacobi-type 
differential equation. As shown below, its normalizable 
solutions are known to be expressed in terms of Jacobi 
Polynomials P^^'^iC] (sec Appendix C for our conven- 
tions). Changing the dependent variables as 

a™.[C] - (1 - C)^l™'^l(l + C)^I"+^I/3™.[C], (44) 

and using the fact that m is half-integral [Eq. (P7|) ]. one 
can verify. 



\m\(\m\ + 1) + - \ 



(45) 



Comparing this with the standard form of the Jacobi's 
differential equation (c./. Eq. (jC3|) ). 



d_ 

dC 



(1 - C')^ + {'^ - - (m + + 2)C} - 



one can identify the parameters as 



(46) 





(T 


TO 




a 




~ 2 


\m\ 


TO — 


2 




cr 


m 




a 


m| 


+ 2 


\m\ 


m + 


2 



(47) 



and 



A^ = n(n + /i + 1/ + 1) + |to|(|to| + 1) + ^ 



n + \m\ + 



(48) 



where the normalizability of the wave function requires 
n to be non-negative integers. 



71 = 0,1,2,3, ••• 



(49) 



In Eqs. (|47| cr = ± referes to the subscript of l3,na[C]- 
The second equality of Eqs. (|47|) holds since m is half- 
integral. Thereby, apart from a normalization constant 
Cnrrun which will bc determined later, the surface wave 
function /3m<T[C] is expressed in terms of the n-th order 
Jacobi polynomial as 



/^mo" [C] — /^nmcr [C] — ^nmcrRn " " [C] ■ 



(50) 



Clearly, n is common to /3m+[C] ^-nd /3m-[C] once A is 
chosen to be an appropriate quantized value imposed by 
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Eq. (gSl). Substituting Eq. (gS]) into (gT]), one finds the 
surface energy spectrum, 



A ( 1 

E^^W in+\m\ + - 



(51) 



where m and n take, respectively, half-integral and non- 
negative integral values [see Eqs. (|57|) and (jH])]. The 
energy spectrum ([5T|) has a couple of specific features: 

1. The discrete energy levels E\ are placed with an 
equal distance, A/R, symmetrically on the positive 
and negative side of i?, excluding the zero energy, 

£; = o. 

2. The degeneracy g\ of each energy level increases 
linearly with \E\, taking every (positive) even num- 
bers, g\ = 2, 4, 6, • • • , when one starts counting it 
at the first positive and negative energy level with 
A = ±1. 

The discrete energy levels E\ can be also casted from 
the view of finite-size energy gap. On the surface of an 
infinitely long cylinder the surface electronic spectrum 
shows an energy gap as a consequence of the spin-to- 
surface locking and the resulting half-integer quantiza- 
tion analogous to Eq. (P7|) . The size of the obtained 
energy gap on a cylindrical surface is inversely propor- 
tional to the radius of the cylinder. Here, in Eq. ([5T|) 
the zero energy E = Q is indeed excluded, which can be 
regarded as a remnant of the energy gapped in the con- 
tinuos spectrum. Note that this has nothing to do with 
the discreteness of the spectrum due to finite-size "quan- 
tization". What counts here is the absolute position of 
the entire (discrete) spectrum. The size of the "energy 
gap", 



A(i?) =Ei- E^i 



2A 

the radius 



(52) 



of the 



is again inversely proportional R, 
spherei^ 

In order to determine the relative magnitude and phase 
of Cnm+ and c„m_ , one needs to go back to Eqs. ((39| and 
substitute /3mcr[C] obtained as in Eqs. ((5(1)) into this cou- 
ple of equations. This is straight-forward, but turns out 
to be rather a tedious work. Leaving explicit demon- 
stration of this to Appendix D, here we refer only to its 
result: 

c„m- = -sign[Am]c„,„+. (53) 
Choosing c„„j+ to be real, and taking the normalization. 



27r / d(cos0)|a„™(0)|2 = 1, 



(54) 



also into account (c./. Eq. (jC6p for the normalization of 
Pn^[C,]), one can give an explicit list of the coefficients 



in which 



v/n!(n-|-2|TO|)! 



2^2l™l+i/2(„+ _ 1/2)!- 



(55) 



In Eq. (j54|) we introduced the notation q:„,„(0), i.e., 
"m(^') = OLnm{S) in accordaucc with Eq. ([50)) . 

As a solution of the effective surface Dirac equation 
(|M)). Eqs. dSS]), dSS]), din]), dM]) specify an exphcit form 
of the surface eigenspinor, <y.{9,(f)) = e'™'^Q:„„i(0), with 
a definite half-integral angular momentum m. It may be 
suggestive to give a few concrete examples of OLnm (^) for 
small values of n and |m|. For simplicity let us restrict 
ourselves to the case of A > (positive energy). For 
n = with an arbitrary half-integral m, a„m(0) is found 
to be 



OLarn{d) = 



(56) 



1 y(2H)! 

i)' 



47r TO 



(sin 1)1"- 1 (cosf) 

/ . 0\\m+^\ I p\ 

sm ^ cos ■ 



|m-i| 



2 / V^"'^ 2 J 

where Sm is an abbreviated notation for 



sign[TOj 



TO 

ItoI 



More specifically. 



Oinl = 



«0i + 




1 


cos 1 






V47r 


_ - sin f _ 



An 
2^ 



sin f cos^ I 
' sin^ I cos I 



For n = 1 and to ~ ±1/2, 



1 



2V27r 
1 



(1 



sin 1(1 



3 cos 9) 
3 cos 6) 



2V27r 



sin 1(1 + 3 cos 6*) 
cos |(1 — 3 COS0) 



(57) 

(58) 
(59) 
(60) 
(61) 

(62) 
(63) 



The two eigenstates ()58p and ([55)) constitute the two low- 
est energy (degenerate) eigenstates on the E > side, 
whereas Eqs. ()60ll63p correspond to the four second low- 
est (first excited) states. 

Let us point out that the surface spinor wave functions 
thus constructed imply a rich spin texture that the sur- 
face eigenstates manifest. This can be already seen in the 
simplest example, the case of the eigenstate cxqi given as 
in Eq. ((58)). At the north pole {6 = 0) the spin in this 
eigenstate is pointed in the -|-2;-direction, perpendicular 
to the tangential plane of the sphere at this point. As 
6 increases, the spin tends to lie closer to the tangen- 
tial plane. At 9 = tt/2, i.e., on the equator the spin is 
completely in-plane to the spherical surface. Now, stay- 
ing at 9 = tt/2, if one let (j) vary from to 27r, i.e., as 
the electron hypothetically travels around the equator. 
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FIG. 1. ^-dependence of the two spin components of a^i{6). 
Q?i i + and aj^i_(S), shown, respectively, in red and in blue. 
The two spin components acquire an equal weight at — 9i ~ 
0.2327r, 9 — tt/2 and at 6 = tt — 6i. a-^^i_^{9) vanishes at 
9 — 92 — 0.3927r and at 9 = n, whereas a-^i_{9) vanishes at 
9^0 and at 9 ^ iv - 62. 

the spin also completes a 27r rotation in the (a;, j/)-plane, 
manifesting the feature of "spin-to-surface locking". As 
one further increases 9, the spin down component starts 
to dominate, before it dominates completely at the south 
pole. 

Another example showing a more complicated spin tex- 
ture is the case of aii{6) given as in Eq. (|62|) [the two 
spin components of aji (0) are depicted in FIG. [Tj. At 
the north pole the spin in this eigenstate is again pointed 
in the -l-z-direction. As 9 increases, it starts to lie, but 
tends to lie more strongly than the case of otQi . At 

9 = 2 arecos ^-y^ - 0.2327r = 9i , (64) 

the two spin components acquire an equal weight. At this 
angle 6* = ^i, if one let <j> vary from to 27r, complete 
spin-to-surface locking occurs. As 9 exceeding 9i, the 
|f— )-component starts dominate. At 

e = arecos[l/3] ~ 0.3927r = 6*2, (65) 

the centrifugal spin component q:ii^{9) vanishes. There- 
fore, at this point the |r— )-component dominates com- 
pletely, and the spin is pointed to the center of the sphere. 
As 9 further increases, the spin gradually tilts back to the 
tangential plane (at 9 = 7r/2, spin-to-surface locking is 
recovered), then it starts to be further tilted toward the 
outside of the sphere. At 6* = arccos[— 1/3] = tt — 02, the 
centripetal spin component a-^i_{9) vanishes, and the 
spin finds itself purely in the \r+) state. At 6* > tt — ^2, 
|r— ) starts to dominate again, and the spin is finally 
pointed in -l-z-direction at the south pole. The behavior 
of the spin as varying (j) from to 27r is the same as the 
case of cxqi since the two states share the same quantum 
number m = 1/2. 



FIG. 2. Pseudo-degenerate spectrum E[J] of the surface 
eigenstates, obtained by diagonalization of the tight-binding 
Hamiltonian on a cubic lattice. The system has a cubic shape 
of linear dimension L = 16. mo/m2 = —1, A = 1. 



The spin rotates more drastically also in the 4>- 
direction for the eigenstates with \m\ > 3/2. Clearly, 
the surface spinor wave functions with quantum numbers 
n and |m| higher than the examples given in Eqs. (j58t - 
I63p show a richer spin structure on the spherical surface 
as a consequence of the interplay between the two types 
of Berry phase; c.f. the case of cylindrical geometry in 
which only a single type Berry phase manifests. 



V. COMPARISON WITH THE 
TIGHT-BINDING CALCULATION 

We have so far investigated specific features of the 
topological insulator surface state occupying a finite vol- 
ume, taking as an example the case of spherical geometry. 
Starting with the gapped bulk effective Hamiltonian, we 
have derived and solved the surface Dirac equation, from 
which we have deduced the surface electronic spectrum 
(|5ip and the explicit form of the spinor wave functions [as 
given through Eqs. p4l[50l[53l[55l) ] on a perfect spherical 
surface. The role of two distinct types of Berry phase has 
been revealed. Here, we take another viewpoint; namely, 
we go back to the bulk effective Hamiltonian [Eq. ([1])], 
and implemented it as a nearest-neighbor tight-binding 
model, which allows for obtaining the spectrum of the 
surface solutions by exact diagonalization. We show that 
the basic features on the surface energy spectrum we have 
found so far in the idealized spherical geometry with ex- 
act rotational symmetry is still valid when that symmetry 
is weakly broken. 

Let us employ the following lattice implementation of 
-ffbuik, i-e., Eq. dD) on a cubic lattice: 

^^lattico = T^m{k) + At^ ^ a J sin kj, (66) 

j=x,y,z 
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where 

m{k)=mQ + 2 ^ 7712(1 — cos /cj), (67) 

is a lattice version of Eq. ([2]). The model specified by 
Eqs. and (p7)) can be regarded as a tight-binding 

model with only the nearest neighbor hopping. This cou- 
ple of equations determine the structure of energy bands 
over the entire Brillouin zone, which also reproduces, in 
the vicinity of the F point, the bulk effective Hamilto- 
nian, -ffbuik [Eq. ([1])] in the continuum limit. The system 
we consider here has a cubic shape of linear dimension 
L, in which the lattice points placed with a unit lattice 
spacing are restricted to 

l<x,y,z<L. (68) 

The obtained surface energy spectrum E ~ E[J] is 
shown in FIG. [2] Here, J is an index for numbering 
each surface cigenstatc with increasing order of E. For 
an aesthetic reason, and also to case counting of the de- 
gree of (pseudo-) degeneracy we have chosen this index 
J to be half integral, J = ±1/2, ±3/2, • • • . Number- 
ing of the eigenstate is done in such a way that starting 
with E = J = 1/2,3/2,5/2, • • • is attributed to each 
eigenstate with increasing (decreasing) order of E on the 
positive- (negative-) E side [i.e., a positive (negative) J 
corresponds to a positive- (negative-) energy level]. De- 
pending on the number of lattice sites contained in the 
system (= L^) and on the value of model parameters, 
certain numbers of states appear in the bulk energy gap. 
Those surface eigenstates of relatively small \E\ < mo 
are also characterized by the spatial profile of the corre- 
sponding wave function; their wave function is localized 
in the vicinity of the cubic surface. In this demonstration, 
the system's size is i = 16, and the model parameters 
are chosen as 7710/7712 = —1, A= 1. 

E = E[J] plotted in FIG. [2] shows a "pseudo- 
degenerate" feature much reminiscent of the quantiza- 
tion characteristic on the spherical surface; the spectrum 
([CT]) . and the degeneracy rule [Eqs. (gHl) and (1^]. No- 
tice that the zero energy state i? = is clearly excluded. 
Horizontal gridlines are located at the positions of E that 
are an integer multiple of £'[1], the first (positive) energy 
level. These gridlines are shown for verifying that the 
energy levels are equally placed, one of the character- 
istic features of the spectrum on the spherical surface. 
Small deviation from the "expected" spectrum ([STj) can 
be seen for |A| > 3, which can be interpreted as a con- 
sequence of the breaking of spherical symmetry. Verti- 
cal gridlines are drawn for highlighting the degree of the 
pseudo-degeneracy of each level. If one recalls that a 
magnetic monopole is effectively induced at the center of 
the sphere, the (pseudo-) equally spaced spectrum illus- 
trated in FIG. [5] can be interpreted as Landau levels. For 
X = 16 the value of E[l] is found to be E[l] ~ 0.0880 in 
units of hA/a, where A is the group velocity of the sur- 
face state, and a is the lattice constant. Taking experi- 
mentally realistic values for A and the characteristic 



energy scale HA/a is estimated to be on the order of 0.1-1 
eV. 



VI. CONCLUDING REMARKS 

The protected surface state of topological insulator has 
an "active" property that it reveals only when it is em- 
bedded onto a curved surface. On a cylindrical surface, 
it induces an effective magnetic solenoid of total flux ±7r 
In the same sense, a magnetic monopole of strength ±27r 
is induced when it is embedded onto a sphere. In this 
paper we have explicitly examined this active property 
of the topological insulator surface state, focusing on a 
most suggestive case of the spherical surface. As a re- 
sult, the following unique profile of such surface states 
has been found. The two important features are 

1. a unique quantization rule; equally spaced spec- 
trum with the exception of the "evaporated" zero 
energy state, and the simple degeneracy rule, 

2. a rich spin texture resulting from the nature of com- 
plicated spinor wave function. 

These characteristic features derived analytically using 
an idealized spherical geometry is then contrasted with a 
tight-binding calculation on a cubic lattice of cubic shape. 
Unexpectedly profound agreement of the two results sug- 
gests that those features which we have demonstrated 
on a perfect sphere capture the essential characteristic 
of the surface states occupying a finite volume of more 
generic shape, inevitably involving a curved surface. In 
this sense it is natural to expect that the obtained spec- 
trum be applicable to the spectroscopy of topological in- 
sulator nano-particles. In particular, we predict a unique 
photo absorption/emission spectrum resulting from the 
equally spaced energy levels of the low-energy surface 
eigenstates. 

Most of the existing works characterizing the topolog- 
ical insulator surface states are based on Bloch states. 
Topological properties, however, do not depend on the 
translational symmetry. Here, we have demonstrated 
this, focusing on the angular momentuni^"— (instead 
of the linear momentum) as the good quantum number. 
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Appendix A: Symmetry class and topological 
invariants 



In the classification of the Dirac Bogohubov-de Gennes 
type Hamiltonian in terms of the time-reversal (Q), 



our 



particle-hole (C) and chiral (Fs) symmetries 
starting bulk effective Hamiltonian ([T|) falls on the class 
of All, to which Z2 topological insulators in two (2D) 
and three spatial (3D) dimensions are classified. This 
class of models has the symmetry, 8^ = — 1, = 
and = (here, "0" indicates that the system does 
not possess that type of symmetry), and are character- 
ized by Z2-type bulk toploogical invariants both in 2D^ 
and 3D.— For the specific choice, e(p) = 0, this sym- 
metry is upgraded to the class Dili, yielding = —1, 
= 1 and F5 = 1, where for the specific Hamiltonian 
^ C and F5 are given by C = (JyTyK and F5 = r^. This 
symmetry class obeys to a Z-type bulk topological clas- 
sification, characterized by a winding number M for the 
Wilson-Dirac-like operator in three dimensions. 

In the following, we explicitly construct and evaluate 
this Z-type winding number A/". Here, to ease the nota- 
tion we rewrite Eq. ([T]) in the specific case of e(p) = 
as 

H = Arifff^pf, + T3m{p), (Al) 
where = 1,2,3 and m{p) = toq + "^2P^- Note that 

H'- = Ay. + m\p)^R\p)l. (A2) 
Therefore, a deformed Hamiltonian 



R 



(A3) 



has eigenvalues ±1. To characterise the topological prop- 
erty of H, let us make a rotation in the r-space such that 
Ti — >■ T2, T2 — >■ Ta, Ta — >■ Ti. Then, the Hamiltonian is 
converted into 



H 



R \ym{p) + lAcF^^p^ 
( QiP) \ 

V QHp) / ■ 



Here, a 2 x 2 SU(2) matrix Q emerges, 



(A4) 



Q 



R 



(A5) 



satisfying Q^Q = 1 and det Q = 1. Since n3{SU{2)) = Z, 
Q should be characterised by an integer winding number 
gencrically. 



The winding number A^ is defined by 

247r2 
1 



M^TT^ I tr {QUqY 

d^P e^-ptr {Q^d^Q){Q^dM){Q^dpQ), 



247r2 



(A6) 



where d stands for the exterior derivative with respect to 
p^, and df^ is the derivative with respect to p^. For the 
normalization of the above equation, see, e.g., Eq. (66) 
of Ref.i^. 

Let us look into the nature of this winding number 
more precisely. Since Q is a function of p^, the winding 
number Af characterizes the mapping from R"^ to SU(2). 
However, at the infinity of M'^, namely, \p\ 00, Q be- 
comes a single element of SU(2), Q sgn(TO2)l- There- 
fore, in the present mapping, the infinity |p| — >■ 00 can 
be regarded as a single point, which make it possible to 
regard M.^ as S^. By the use of the fact that SU(2)- S^^ 
the winding number A^ characterises the mapping from 
S'^ to S'^, which can be classified by 7r3(§'^) = Z. 

It is not difficult to guess the winding number in the 
following way. At the origin \p\ = 0, Q = sgn(mo)l. 
Therefore, if sgn(mom2) = 1, Q can be defomed into 
Q = sgn(TOo)l by taking the limit |too| ~^ 00, giving rise 
to a trivial winding number M ~ Q. On the other hand, 
if sgn(mo'7i2) = —1, Q at \p\ ~ cannot be deformed 
into Q at \p\ = 00, so that this case gives N = ±1. 

Let us compute the winding number M concretely. Af- 
ter a tedious but straightforward calucation, one finds 



12 



m{p) 



We reach, therefore. 



-12 



(Fp- 



i?4 



Too - m2P 



(A7) 



/R3 [P^ + ("T-0 + TO2p2)2]2 

-12 X (At:) f°° p^{mo - m2P^) 



247r2 



dp 



[p2 + (mo -I- TO2p2) 



2\2]2 



[sgn(mo) - sgn(m2)] 



(A8) 



The last formula indicates that the system is indeed in 
the topologically non-trivial {J\f = ±1) phase when toq 
and 7772 have the opposite sign, albeit in the trivial {M = 
0) phase when they have the same sign. 



Appendix B: The zero-energy condition 

The radial eigenvalue problem considered in Sec. II 
has two basic ingredients: the eigenvalue equation ([5]) 
and the boundary condition ([7]) at r = i? (on the sur- 
face of the sphere). Here, we prove that the solutions of 
this radial boundary problem satisfies automatically the 
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energy condition, E± = 0. This observation paves the 
way for constructing the basis eigenspinors given in Eqs. 
(fT7|) . In other "simpler" geometries, such as the case a 
semi-infinite system with a flat boundary;^ or a cylin- 
drical system infinitely long in the axial direction;^ the 
scenario applies, but the explicit use of the zero-energy 
condition may not be indispensable because of the sim- 
plicity of the problem. 

We first need to find the explicit matrix form of H± = 
^^buik|pe=o,p^=o- Focus On the three momentum opera- 
tors, p± = —id±, Pz ~ —idzi that have appeared in Eq. 
([3]), and in addition, = — in m{p). In the spherical 
coordinates, these can be expressed as 



d± 



■ n9 1 

sm6l— -I- - 
or r 

e±'^ sme—+D±, 
or 



cos 



89 sine 



d 



(Bl) 



d - cos 6— - ^ 
^ dr r 89 

d 

= cos 6*— + 
or 



and 



dr 



2d_ 

r dr 



I? 



(B2) 



(B3) 



where is a square of the orbital angular momentum 
operator. Since D±, and involve only angular 
derivatives, we put them in H\\ . Keeping only the first 
terms of Eqs. JET]), (|B2| and ((B3| . one finds. 



—ina 


—inbe^'^ 









m± [k] 
iKa 



-inbe-"'' 




(B4) 

where we have assumed a surface solution of the form of 
Eq. ([6]), and replaced the r-derivatives with k, the in- 
verse of the penetration depth. Eq. ()B4|) can be indeed 
regarded as a 4 x 4 c-number matrix specified by a pa- 
rameter K. We used the notation, H±[k], to make this 
point explicit. We have also introduced, for shortening 
the expression, the notations, a = Acos9, b = ^sin6' 
and 



mj^[K] = mo — ?7!,2 [ K H K 



(B5) 



Here, in Ec^. (|B5P an r-depcndence in the last term looks 
cumbersome. But as far as the surface wave function is 
well localized in the vicinity of the surface at r = i?, one 
can safely replace this coordinate r by a constant R. On 
the other hand, as far as the same assumption is applied 
this last term itself becomes negligible, since as far as 



K 



term. 



<C R; the second term is much larger than the last 



Thanks to the symmetric structure of the matrix form 
of Eq. (|B4p . the secular equation, det — E±\ = 

for the radial eigenvalue problem ([5]) becomes as simple 
as 

dct\Hj_[K]-Ej_\ = {K^A^-m^[Kf + El}^ =0. (B6) 

This can be regarded as a quadratic equation for under 
the approximation of m^[K] = toq — m2K^; with its two 
solutions Kj_ satisfying. 



2 .2 _ 

K ( ^ — 



Ei 



(B7) 



Now, in order to cope with the boundary condition ([7|), 
two surface solutions of the form of Eq. (O, one with 
K = and the other with k = k- must be superposed, 



IV') 



(B8) 



where u[k\ is an eigenvector of i?_L[K] given in Eq. (jB4|. 
The only way that this solution be compatible with the 
boundary condition ([T]) is to have simultaneously c+ -I- 
c_ = and = i.e., takes the following 

form. 



(B9) 



where is a normalization constant. 

In the reduction from Eq. (|B8p to Eq. (|B9[) . the sec- 
ond condition stating that two eigenvectors belonging to 
different k's should coincide, was crucial. Indeed, this 
coincidence occurs only under a very specific condition. 
In order to clarify this point, let us consider the following 
quantity. 



_ H^[k]-E^ 



K 

Mi[k\ —ia 

—ia —M2[k\ 

-ibe"'' 

-ibe^'t' 





MM 
ia 



-ibe-^'t' 



ia 
-M2H 



(BIO) 



where 



MM 



m_L [k] — E± 



M2M 



mj_ [k] + El 



(Bll) 



Notice that u[k] is a zero-eigenvalue eigenvector of this 
matrix, i.e., A[k]m[k] = 0. In order that two of such 
eigenvectors u[k\ belonging to different k's (= n±), and 
therefore, to different A[k]'s (A[k_|-] and A[k_]) coincide, 
both A/i[k] and Af2[K] must coincide. Namely, one must 
have, simultaneously, Afi[K_|-] = A/i[k_] and il/2[K+] 
A/2[k_]. Clearly, a solution of the form of Eq. (|B9|) is 
meaningful only when k,+ 7^ K-. Therefore, A/i[k_|_] = 
A'/i[k_] signifies. 



K+K- 



mo 



■E, 



m2 



(B12) 
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whereas M2\k+\ = M2[k_] leads to 

mo - E± 



Appendix D: Proof of Eq. ((53|) 



m2 



(B13) 



Recalling that Eqs. (|B12|) and (|B13p must follow in- 
dependently, one can convince oneself that the surface 
solution must satisfy the zero-energy condition, 



= 0. 



(B14) 



Note that Eqs. (|B12p and (|B13|) are consistent with Eq. 
(|B7p . but impose a stronger constraint on the values of 
k's and i?_L. 



Appendix C: A brief reminder on the Jacobi's 
polynomials/differential equation 

• The explicit form of the Jacobi's polynomials is 
given by the following (differential) Rodrigues' for- 
mula: 



where 



(-1)" 1 d" 



2"n! p^u[C\dx 

Pp.[C] = (i-cr(i-c) 



[(i-C')>M.[C]], (CI) 



(C2) 



As is clear from its construction, Eq. (|C1|) can be 
also expressed in the form of a contour integral (the 
miegra/ Rodrigues' formula). 

The Jacobi's differential equation is a simple 
rewriting of the hypergeometric differential equa- 
tion. 



e(i-o^ + {c-(i+^+^)e}| 



-AB 



F(yl,S,C,e) = 0, (C3) 



by the change of the independent variable, 

^ 2 ' 
and choice of the parameters, 

A ~ —n, 

B = n + ^ + h'+l, 
C = M+1. 



(C4) 



(C5) 



• The Jacobi's polynomials satisfy the follow- 

ing ortho-normal relation. 



rfCPM^[C]PC[C]^C[C] (C6) 

2^+'^+ir(ni +fi+ l)r(7ii + 
ni!(2ni +fj, + v + l)T{ni + v + 



In order to determine the relative magnitude and phase 
of Cnm+ a-nd c„„i_, one needs to go back to Eqs. (|39p and 
substitute /3,„o-[C] given in Eqs. ([50]) into this couple of 
equations [naturally, the change of the dependent vari- 
ables must be taken into account; see Eq. (HH)]. Chang- 
ing the independent variable from to = cos 0, let us 
rewrite Eqs. ([5^ as 





m + 


C/2 


dc 


Vl 






m — 


c/2 


dC 


Vl 





a„+[C] =Aa,„_[C].(Dl) 



Performing the derivatives explicitly, and changing the 
variables from a,na to /3,„(j, using Eq. (j44p . one finds 



d 
'dC 
d_ 

'dC 

for m > 1/2, and 

d 
'dC 
d_ 

'dC 



(1-0- 
(1 + 0- 



(1 + 0; 
(1-0- 





















-0] 







A/3„_[C] (D2) 
A/3„_[C] (D3) 



for TO < 1/2. Recall that pma^s as given in Eq. ((50)) are 
proportional to the n-th order Jacobi's polynomial. Eqs. 
(|D2[) and (|D3[) can be further simplified on account of 
the following identities [c.J. Eqs. (A. 7a) and (A. 7b) of 
Ref.— ] , applicable to the derivative of the Jacobi's poly- 
nomials with a specific choice of parameters /i and that 
are implied in these relations through Eq. (|50p . i.e.. 



dC 
dC 



(D4) 



n + m+ - ] Pn 



,Tn+i 



(D5) 



1 



n + TO + - 1 P„ 



where to > 1/2. These identities can be explicitly ver- 
ified, e.g., by the use of the integral counterpart of Eq. 

m- 

The final part of the proof of Eq. ((53)) lies in the com- 
parison of Eqs. ((D2|, ((Ell and ((Dl]), ([D5]). For to > 1/2, 
one can safely take off the operation of absolute value 
to the superscripts of Jacobi's polynomial in Eq. ()50p . 
yielding 



j^vna [C] — l^nma [C] 



P, 



"[C]. (D6) 
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Then, by simply comparing Eqs. (|D2[) with (jD4[) and 
(|D5p . and recalhng A = ±(n + m + 1/2), one can verify 
|cnm+| = |cmn-| with a rclativc sign of 



Cnm- = -sign[A]c„m+. 
For m < —1/2, notice that 



a 

m 

2 



i.e., for such m Eq. ([SH]) becomes 



/^ma"[C] — /^nmo" [C] ^nraaPi 



|Tn| + f ,|m|- 



"[C] 



(D7) 



(D8) 



(D9) 



This allows for the use of Eqs. (jD4l) and (jD5|) with m 
replaced by \m\ in the couple of Eqs. (|D3p . Taking note 
of A = ±(ri + |m| + 1/2), one can again verify |c„„i-|-| = 
|c„„i_|, but this time with a relative sign of opposite 
value. 



sign [A] Q 



(DIO) 



The relations (jD7[) and (|D10p . respectively, for the two 
possible regimes of m complete the proof of Eq. ([55]) . 
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